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I. INTRODUCTION 


We observed Jupiter with the Westerbork Radio Telescope at a frequency 
of 1412 MHz in December 1977 when it was at its most northerly opposition 
to the earth. As was pointed out in Ch. Ill of the previous final 
technical report for grant NSG 7264 (hence forth referred to as Paper I) 
that was an optimal time to observe this planet. 

We obtained pictures of Jupiter in all four Stokes parameters at 24 
different rotational aspects of the planet, each integrated over 15 degrees 
of Jovian rotation. The maps clearly indicate the presence of higher order 
terms in the dipolar field of Jupiter at distances of ^ 2 Rj from the 
center. They also show a displacement of the main dipole of 0.119 + 0.009 Rj 
from the center of the disk towards longitude 135°-145° (Syst. Ill 1965.0) 
and a displacement of 0.04 ± 0.04 Rj towards the north in agreement with 
the displacement found by the Pioneer spacecrafts . From the data I would 
estimate the thermal disk temperature at this frequency to be more than 
300 K but less than 340 K, which implies an ammonia mixing ratio of m 5x10 
(Gulkis, 1973). 

II. OBSERVING AND REDUCTION TECHNIQUE 

The properties of the Westerbork telescope and its operation have been 
described extensively in Paper I. Since the circularly polarized flux of 
Jupiter can give much information concerning the magnetic field of the 
planet we tried to map this flux with an accuracy much better than generally 
obtained. Under normal circumstances the dipoles of the Western telescopes 
of the interferometers are aligned at position angle PA = 90° and those 
of the Eastern telescopes at 45° so they are "crossed" to each other (see 
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Fig. !). If the dipoles are set parallel to each other it appears that we 
can map the circular polarized flux V with a much better accuracy, together 
with either U or Q depending on the position angles of the dipoles. 
Following the papers by Weiler (1973) and Weiler and Raimond (1976) we 
see that the complex outputs R of. the four channels of an interferometer 
with the "crossed" configuration can be written as: 


= j pi 1 c «(' < 1 - -Q{i+d: x - iQ- x ) 

+ 1/ { i — a xx + .•«;,} — » v { 1 + a; x - /e; x }) 

R » = 2]72 ' G *yd U+4T,- iO - e { 1 - 4;, + } 

- V { 1 + j;, - «0; y } + tv { I - j; y + }) ^ 

= 2 ^k Gyx(I { 1 + dyx ~ ,6 >* } + <20 - K + m;,} 

+ U{i+a; x - io; x ) + « v { 1 - a; x + io; x }) 

Ryy = J\n Gyy(l { 1 “ Ayy + i6yy} + <2{l+^yy- i0; y } 

- U { 1 - a;, + id;,) - iv { 1 + A - - io; y ) y 


where the G's are the complex instrumental gain factors, I, Q, U and V - 
are the complex visibility functions of the four Stokes parameters. 

A = Ay - Ap are errors in the dipole setting of the Eastern and Western 
telescopes, i.e. an instrumental linear polarization, and 0 = 0y + 0^ 

represents an instrumental circular polarization. Both terms A and 0 can 
be adjusted to a level of < 1% of I. Neglecting these small terms in the 
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Fig. 1 Dipoles for the "crossed" resp. "parallel" configuration (from 
Weiler, 1973 and Weiler and Raimond, 1976). Each telescope 
has the two dipoles x and y. 


solution for the four Stokes parameters cause errors of the order of 
0.4% of I in the individual maps. Fluctuations in the system gains are of 
the order of £ I % of the receiver output giving a total error of ^ 0.5 ^ 
of I in all maps. 

Setting the dipoles of all telescopes parallel to eachother at position 
angles 90° gives us the equations: 


* iQ- u\a; k - /o M ) - /f(/i vv i 

• -><0+Q{A.: r -/t) xy )-U +,v\ 


I 

i 
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N » = iGrAl + Q+UM^-so-i-sVM^-nty) 
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Neglecting the terms (U, Q, V) (A, 0) in this configuration cause errors 
of 0.1 % of I in the various maps and we can write the receiver output as: 


CM 

2<? a , l.l I / |/| 

/ 

l-V). 
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where e = (A - 0 + ), a complex instrumental error which can be found from 

calibration measurements as will be shown below. Together with the fact 

that possible gain fluctuations in the system are < 1% of the receiver 

signals it will be clear that the disturbances in U and V generated by 

system errors causing "leakage" of I will be 'v 0.1 % of I or less, 12 hr 

observations of the planet were made on each of five consecutive days in 

order to get finally two dimensional pictures of Jupiter at all rotational 

aspects, smeared over only 15° of central meridian longitude (Dec. 16 - 21, 

hr 

1977; one extra observation of 7 l was made on the first day). To be able 
to calibrate the data we chose an observing sequence as listed in Table I. 
with the calibration sources shown in Table II. 

After all measurements were corrected for known instrumental and sky effects 

(Paper I and van Someren Greve, 1974) we used the adopted flux of 3C147 

to define the gain of the system assuming zero polarization for this source. 

The other three sources listed were used as geometry and phase calibrators. 

The complex correction factors G and G could be determined directly 

xx yy 

from the calibrations taken in the "parallel" configuration. However, we 
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Table 1 : Observing sequence 


Sidereal day 

350 

351 

352 

353 

354 

355 

356 

357 


Source 

3C48 

3C286 

3047 

Jupiter 1 

3C286 

3C48 

3C147 

Jupiter 2 

3C286 

3C309 . 1 

3C48 

3C147 

Jupiter 3 

3C286 

3C309 . 1 

3C48 

3C147 

Jupiter 4 

3C286 

3C48 

3CI47 

Jupiter 5 

3C286 

3C48 

3047 

Jupiter 6 

3C286 

3C286 

3C48 

3047 

3C48 

3047 

3C48 

3C286 



Table II : Calibration sources. The flux densities of 3C48, 3C286 
and 3C309.I are "mean" flux densities, derived from 
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are left with two unknowns in both the r' and r' . Using the calibration 

xy yx 

observations taken with the "crossed" configuration we were able to 

define the G and G for each observation due to the fact that the 
xy yx 

mutual channels of one interferometer are very stable (fluctuations are 

less than 1% in amplitude, and 0?5 in phase). The observations of the 

unpolarized source 3C147 provided the complex instrumental error terms 

£ and £ . The Fourier components for each interferometer of the four 

xy yx 

Stokes parameters can then be represented by: 


I = \ (xx + yy) Q = \ (yy - xx) ORIGINAL PAGE IS 

OF POOR QUALITY 


U = |{-2xy + 2yx + (e - e )xx + (e - e )yy} 
* J J xy yx xy yx y 3 


V = “ 4 ( 2 xy + 2 yx - (e xy + 


e )xx - (e 
yx xy 


+ e )yy } 
yx 


where xx, xy, yx and yy are the calibrated complex receiver outputs. It 
appeared that the maximum circular polarized intensities detected in the 
measurements of 3C286 and 3C48 was less than 0.1% of the maximum intensities 
of I; thus the instrumental error in the V maps indeed is less than 0.1% 
of I. The degree of linear polarization of 3C286 is 9.4 ± 0.05 % at an 
angle of 117°. These angles are not corrected for Faraday rotation in the earth' s 
ionosphere which will cause errors of less than 'u 5° (Weiler and Raimond, 1976). The 
polarization of 3C286 is in excellent agreement with the numbers published by Berge 
and Seielstad ( 1 972) (P = 9.36 ± 0.07 , PA = 32 + 1.1) and Morris and Berge 
(1964) (P = 9.3 ± 9.5; PA = 32 + 2); the first publication yield numbers 

L 

for 3C48 of -P = 0.43 ± 0.07, PA = 139.7 ± 4.8, the second one of 

P^ ?= 0.6 + 0.6, 0.7 + 0.6 and PA = 154 + 28, 162 ± 25. Morris and Berge 

(1964) also measured 3C147 to have a P of 0.6 ± 0.3 and 0.2 + 0.5 with 

L 
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PA = 173 * 10, 138 t 57. The position angles are not corrected 

for Faraday rotation in the earth's conosphere. Since observations of 
3CI47 and 3C48 show very small values for the degree of linear polari- 
zation it will be clear that these numbers and in particular the posi- 
tion angles largely depend on the "chosen" polarization of the cali- 
bration source. 

At its time of observation Jupiter was at a distance of 4; 1525 AU 

from the earth, so that the measured flux densities were multiplied by 
4 1525 ^ 

a factor of ( ~ - ) to correct them to those at the standard distance 
4 . 04 

of 4.04 AU. 

After removing the effect of the planetary motion along the sky from 
the Jovian data (as we did in Paper I) , we cut the observations into 
25 min. pieces, which is about 15° of Jovian longitude. This creates 
a smearing effect on the data which amounts to ^ 20 arcsec at an an- 
gular distance of ^ 3 Rj . 

Fig. 2 shows very schematically the coverage of the UV plane 
for these data together with a picture of the beam with which we have 
"observed" the planet. Compared to the 1973 data presented in Paper I, 
we now have a much better coverage of the UV plane (6 observations 
equally spread in hour angle) and much less trouble with the quality 
of the calibration since the planet was at a much higher declination 
and thus elevation. 

Using the procedure called CLEAN, the deconvolution technique 
developed by Hogbom (1974) we were able to remove the antenna pattern 
from the maps so that after having restored the components to the map 
using a Gaussian shaped beam (with a half power width of 1 6 6 ) the un- 
certainties in the maps appeared to be ^ 17. of the maximum intensity 
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Fig. 2 The upper part of the figure shows very schematically the 

filling of the Fourier transform (UV) plane used to map Jupi- 
ter at all rotational aspects of the planet. The lower part 
shows the beam with which we are looking at the planet at a 
particular longitude. 


2 2 

in the I maps, 1.5 - 2.0% of the maximum response in the polarized (/Q +U ) 
maps, and the V maps were accurate to 0.1% of the maximum intensities of the 
I maps. The position angles were corrected for Faraday rotation in the iono- 
sphere using electron density data obtained from the Dutch meteorological 
Institute in the Bilt, which gave a correction of 1.43 with variations up 

Or 

to 20 - 30%. The error remaining in the position angles is 0.5. 
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III. INTEGRATED JOVIAN PARAMETERS 


Fig. 3 and 4 show various parameters of the integrated Jovian radiation 
as functions of, respectively, Jovian magnetic longitude (Syst. Ill 1965.0) 
and magnetic latitude of the earth with respect to Jupiter. The parameters 
are calculated by integrating the two-dimensional Jovian maps in I, Q, U 
and V, which all are averaged over 15° intervals. 

Indicated are: 

a_ the total integrated flux density in Jy 

b polarization direction: the position angle of the electric vector measured 
eastward from the north in the sky (corrected for Faraday rotation in the 
earth's ionosphere, see Section II) 

£ degree of linear polarization 
<1 degree of circular polarization 

£ (in Fig. 3) the magnetic latitude of the earth with respect to Jupiter. 

This latitude was calculated according to 0 = D„ + 8 cos (1 XXT - 1 ) 

& m E III o 

(as in Paper I). D^, is the declination of the earth relative to Jupiter's 
rotational equator (+ 2?24 at this epoch), 8 is the angle between Jupiter's 
magnetic and rotational axes taken as 10° and 1 is the central meridian 
longitude of the magnetic north pole taken as 201° (as in Paper I). 

The variations with longitude are consistent with the beaming found by others 
from the integrated data (see Paper I). Since the present observations were 
made in an epoch when the declination of the earth D„ was positive we find 

Cj 

a much deeper minimum at 1 ^ 200° than at 1 ^ 20°. Again the maximum in S 
at 'v 300 is higher than at a, 100°. With the help of the data shown in this 
paper we can interpret this difference as being caused by the presence of a 
region of enhanced emissivity near 1 250° (as was already suggested by 

de Pater and Dames, 1978) together with the fact that the magnetic field has 
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Fig. 3 Various parameters of the integrated radiation from Jupiter 
at 2 1 cm as functions of the magnetic longitude of Jupiter 
(in Syst III 1965.0) 

a) integrated flux density in Jy 

b) position angle of the electric vector measured earthward 
from north in the sky 

c) degree of linear polarization 

d) degree of circular polarization 

e) magnet ic latitude of the earth with respect to Jupiter. 





Fig. 4 The parameters a-d from Fig. 2 are plotted as functions of 
the magnetic latitude of the earth with respect to Jupiter. 







a multipole character rather than a dipolar one at 'v 2 Rj from the. center 
(see section IV). The degree of circular polarization in Fig. 2 clearly shows 
maxima at the longitudes of the magnetic poles with a slightly higher 
maximum at 1 ^ 200° caused by the positive declination D..,. The curve how- 
ever shows a rather large deviation from a pure sine curve, which is 
pronounced in Fig. 4 by the flattening of the curve at the ends of This 
flattening is determined by the pitch angle distribution of the electrons 
(Roberts and Komosaroff, 1976), which in general can be written as I sin^ 
(* e , where a g is the equatorial pitch angle. When q = 0 the distribution is 
isotropic, larger values of q specify flatter pitch angle distributions. 

Roberts and Komosaroff (1976) point out that to fit the beaming of the radiat- 
ion and the curves for linear and circular polarization for their data they 

need a q. = 1 and q„ > 30. 

1 l 2 ^ 


IV. TWO-DIMENSIONAL JOVIAN MAPS 


Fig. 5 shows the two-dimensional maps of Jupiter at 21 cm at all 

rotational aspects of the planet averaged over 15° of Jovian central meridian 

longitude. Indicated are maps of the total intensity, circular and linear 

polarized flux and the magnetic field of the planet. This last map is made by 

rotating plots of the position angle of the electric vector (corrected for 

Faraday rotation in the earth's ionosphere as mentioned in section II) over 

an angle of 90°. It is assumed herewith that there is no Faraday rotation 

within Jupiter's magnetic field, a generally accepted idea which was confirmed 

in Paper I considering the integrated parameters of the source of A = 21 cm 

and X = 50 cm. For individual parts of the source we can write the rotation 

measure R as: R = 2.62 x 10 / n .H.dl where R is in radians/ 2 , n 

mm e — m cm e 

3 

is the number of thermal electrons/^3 (50-100/cm , Intrilligator and Wolfe, 
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1976 - but they point out that it is very difficult to measure this amount!) 

H the magnetic field strength along the line of sight in Gauss (exactly zero 
for a pure dipole field if you are in the equator!) and 1 the distance in 
cm. We see that the Faraday rotation R^A 2 near 1-2 Rj from the center at 
A = 21 cm is less than ^ 2°. The total intensity maps show the existence of a 
region of somewhat increased emission which we shall call a "hot" region, 
imbedded in the Jovian field at a longitude of 'v 250°. However, at some 
longitudes where we expect to see this region at one side of the disk there is 
practically no difference in the amount of flux received from both sides of 
the disk - see f.e. the maps at ljj. ^ 110°; 125°. This same feature is- shown 
in the maps in paper I. 

The maps of the circular polarized flux are not at all symmetric at these 

longitudes. In this report I assume that most electrons are near the magnetic 

equator with pitch angles near 90°. In the future I will find out whether 

it is possible to get changes in this distribution at places where the field 

is distorted and whether or not a change will influence our model. The 

circularly polarized flux is a measure of the magnetic field component 

directed along the line of sight. Its intensity is related to the square root 

of the strength of the field along this direction and its sign shows whether 

the field is directed towards us (left handed polarization 'u negative sign) 

or away from us (right handed polarization 'v positive sign). We therefore 

tried to roughly interpret Jupiter's magnetic field with help of these maps 

of the circular polarization and the total intensity. In Fig. 6 I plotted 

the maximum amount of circularly polarized flux received within a square of 
o 

158 arcsec~ for the left side of the belt as a function of magnetic longitude 
(6a) together with the maximum intensity from that side (6b). Taking a model 
for the radiation belt as shown at the right side with a region of enhanced 
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emissivity of 11-12% at 'v 250° with a sign of 'v 80° in longitudes centered 
at 1 ^ 250° of ^ 11-12% we can predict the variation with longitude 

of both kinds of observed flux density. The sinusoid expected for the 
circular flux is sketched with a thin line in Fig. 6a; the dotted curve in 
b shows the variation in flux expected from only a rotating "hot region" 
in the field; modifying this curve with the sinusoidal variation as expected 
from the beaming caused by the wobbling of the equator gives us the thin 
curve drawn in b. Large deviations between the observed and expected curves 


Fig. 5 Two dimensional maps of Jupiter: indicated are from left to 
right maps of the total intensity, circular and linear po- 
larized flux and a vector diagram of the magnetic field of 
Jupiter (see text). The magnetic longitude of the central 
meridian (Syst. Ill 1965.0) is indicated in the upper left 
corner of each serie of 4. The contour values in mJy/beam are: 
I : 7.5 (lower contour); 50-950 in steps of 100 
V : 1.50-21.75 in steps of 2.25 
P^: 7.5 (lower contour), 25-350 in steps of 50. 

The beam area is 18.093 points. The grid point spacing is 
4 1 8 6 (= 0.177 R j ) in RA and 1 0V636 in Dec. The magnetic 
equator and planetary disk. are indicated in the pictures. 




0° 4Q° 80° 120° 160° 200° 240° 280° 320° 360° » X ffl 11965.01 


Fig. 6 In the upper part the mean maximum amount of circular pola- 
rized flux received within a square of 3x3 gridpoints at the 
leftside of the belt is plotted as a function of magnetic 
longitude (thick line), together with a sinusoid expected 
for a simple rotating dipole field (thin line). The lower 
part shows the maximum intensity at that side of the belt 
(thick Line) together with the expected curve based on the 
model shown at the right side (for more details the reader 
is referred to the text). 

The numbers indicated on the right side are in units of 


5 mJy/beam. 
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show up, e.g., Che distance of 240° rather than 180° between the 
longitudes of the zero values of the circularly polarized flux and 
the rather sharp increase in V and decrease in I near 1 ^ 100°-160°. 
The fact that the circularly polarized flux increases much faster 
than expected indicates that, at that location, the magnetic equator 
is directed downwards with respect to the line of sight as if we are 
looking at the longitude of the north magnetic pole. Modifying our 
sinusoidal beaming curve in I in that manner will diminish the dis- 
crepancy between both curves drawn there. 


CONCLUSION 


ORIGINAL page is 
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Both plots clearly show that the magnetic equator whobbles in a 
non-dipolar way and we have to compare our maps with multipole models. 

In a first approximation I took the model of Smith et al. (1976) based 
on Pioneer 11 data to compare with our data. This model consists of a 
series of spherical harmonics; the quadrupole and octupole terms cause 
deviations from a pure dipole model with the largest anomalies near 
longitudes 50° + n.90°. 

Taking into account a region of enhanced emissivity at X 'U 200°-300°, 
(where the gradient in the magnetic field is rather small) Smith's 
model gives rather good qualitative agreement between the variation 
with longitude of the circular polarized flux and the difference in this 
flux between both sides of the belt. The appropriateness of such a model 
to match the radio data was first suggested by E. Gerard (1976). 

In the next one or two years I will develop more detailed multipole 
magnetic field models and I hope to finally get a three-dimensional 
model of the field coupled with an electron and pitch angle distribution 
in that field. Knowing the field structure and particle distribution one 




Fig. 7 The distance of the centre of the dipole field of Jupiter 

from the centrum of the optical disk are shown as a function 
of magnetic longitude. This distance is measured along the 
magnetic equator, positive towards the right (decreasing 
R. A. ) . 


can define the thermal disk temperature with much better accuracy by 
the two different methods described in paper I. 

Estimating this temperature from the maps shown in Fig. 5 with a "mean" 
dipolar field model for Jupiter, we find a temperature more than ^300 K 
but less than ^340 K . This implies a mixing ratio for ammonia (see Paper 
I) of ^5 x 10 ^ . 

Fig. 7 shows a graph of the center of the dipole field of Jupiter 
as a function of the magnetic longitude. Since the higher order terms 
in the magnetic fieLd decrease much faster going away from Jupiter 
(Quadrup ^ ; Gctupole ^r ~*) than the (much larger) dipolar term we 

can approximate the field at 'v 4-5 Rj as being dipolar. Taking the mean 
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position of this dipole field measured along the magnetic equator we 
find a deviation from the optical center of the planet as shown in 
Fig. 7. It appears that the dipole is displaced by 2V8 t 0'.'2 or 0.119 i 
0.009 Rj towards a longitude of n, 135°- 145°, which is in agreement with 
the displacement found by the Pioneer spacecrafts. A displacement 
of the dipole towards the north of 0.04 - 0.04 Rj can be found from 
the maps, also in agreement with the Pioneer results. 


V. SUMMARY 

A preliminary interpretation of the Jovian observations from Dec. 
1977 shows us: 

a) the presence of higher order terms in the dipolar field of Jupiter 
at distances of ^2 Rj from the center 

b) confirmation of the hot region 

c) a displacement of the main dipole from the center of the planet of 
0.119 - 0.009 Rj towards 1 'v, 135°-145° and 0.04 t 0.04 Rj towards 
the north, in agreement with the displacements the Pioneers found. 

d) a thermal disk temperature at 21 cm of ^ 300-340 K, implying a mixing 

-4 

ratio of ammonia in the Jovian atmosphere of ^ 5 x 10 
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